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Using the Green's function of the 3D heat equation, we develop an analytical account of the 
thermal behaviour of superconducting films subjected to electrical currents larger than their critical 
current in the absence of an applied magnetic field. Our model assumes homogeneity of films 
and current density, and besides thermal coefhcients employs parameters obtained by fitting to 
experimental electrical field - current density characteristics at constant bath temperature. We 
derive both a tractable dynamic equation for the real temperature of the film up to the supercritical 
current density J* (the lowest current density inducing transition to the normal state), and a thermal 
stability criterion that allows prediction of J*. For two typical YBCO films, J* predictions agree 
with observations to within 5%. These findings strongly support the hypothesis that a current- 
induced thermal instability is generally the origin of the breakdown of superconductivity under high 
electrical current densities, at least at temperatures not too far from Tc- 
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I. INTRODUCTION 



In experiments in which superconducting films are 
placed in an environment ("bath") that is thermostat- 
ted at a temperature below their critical tempera- 
ture Tc, and are then subjected to a gradually increasing 
current density J, the electric field becomes measurable 
at the critical current density Jc{Th). It then increases 
ever more rapidly with J until, at the "supercritical" or 
"quench" current density J*(Tb), it jumps to values cor- 
responding to the nonsuperconducting state, the discon- 
tinuity at J* becoming increasingly abrupt as the bath 
temperature is lowered. This phenomenon, which is of 
interest not only for the theory of electrical transport 
in superconductors but also in relation to some of their 
most important applications, is still poorly understood. 
The main mechanisms proposed so far may be crudely 
classified in two classes. One, comprising what may 
be termed current-driven mechanisms, basically invokes 
electrodynamic effects dependent on the microstructure 
of the sample.'^l^ The other invokes heat-driven mecha- 
nisms that are essentially artifactual, postulating that a 
small increase in temperature due to the finite duration 
of measurements triggers a thermal runaway.'^''^ (Further 
discussion of both approaches is available'^') 

A weakness of studies exploring the heat-driven ac- 
count has hitherto been their reliance on results that 
were obtained by numerical methods, the limited scope 
of which somewhat obscures their theoretical interpreta- 
tion. In this paper we address this weakness by develop- 
ing an analytical theory of the thermal stability of high- 
Tc films that explains previous experimental and simula- 
tional results.!^ The theory presented is a full 3D model 
that provides a dynamic equation for the temperature of 
the film as a function of time, together with a thermal 
stability criterion with clear-cut predictions for the de- 
pendence of J* on bath temperature and film geometry. 



Its parameters are those of a homogeneous film material; 
no appeal is made to hard-to-quantify microstructural 
defects, which in some other models play the role of free 
parameters that facilitate good fit to experimental re- 
sults. 

We know of no previous studies that significantly over- 
lap with this work. In particular, the monumental review 
by Gurevich and Mints^*' deals only briefly with the ther- 
mal stability of homogeneous superconductors, and then 
only for thin wires and at the hard superconductivity 
limit (Jc ^ J — Jc), conditions that are far removed 
from those considered here. 



II. BACKGROUND THEORY: HEATING AN 
INFINITE MEDIUM 

Consider a point source embedded at tq — (<^, 77, C) 
in an infinite homogeneous medium that at time t = Q 
has zero temperature. The evolution of the temperature 
field T following delivery of a heat pulse at time t = to 
is governed by the heat equation 
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where D is the thermal diffusivity of the medium. The 
solution i^ 
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For a general heating rate density Q (r, t) , substitution 
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of Fourier's law in the heat balance equation 
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(where q is the heat flux and c the specific heat at con- 
stant pressure per unit volume) affords 



1 f)T 1 . 



(4) 



where the thermal conductivity k = cD. The Green's 
function for solution of Eq. Q is Tq [Eq. 

T{r,t)^ f r^^f^TG(r,t|ro,to)diodVo+Tfc (5) 
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where Tf, is the initial temperature of the medium and 
V is a region containing all (7\),to) for which Q (fo,to) is 
nonzero. 



III. THE MODEL 



A. Constructing the model 



We model the film as an homogeneous parallelepiped 
Vi/2 of dimensions ixwxd embedded in the face of a semi- 
infinite region V00/2 of the same material with no discon- 
tinuity between the film and this substrate (Fig.jl]); and 
we assume that heat generated in the film fiows only into 
the substrate, not into the overlying refrigerant, so that 
heat flux through the face of V00/2 and the free face of V1/2 
is identically zero. This allows application of the method 

V0C/2 are the mirror images of 
y2 in the plane of their free faces (Fig. [21), we 
need only perform calculations for V, the union of Vi 



of image if V1/2 and 



V1/2 and Vo 



and V1/2', and since V is a region of a homogeneous in- 
finite medium, and contains all the sources heating this 
medium, we can use Eq. ([S]) to solve our problem. 

Various features of the proposed model invite justifica- 
tion. In the first place, the assumption that heat gener- 
ated in the film flows only into the substrate is an accept- 
able approximation if the overlying refrigerant is gaseous, 
and also if liquid nitrogen is used, since the thermal con- 
ductivity of the solid substrate can easily be two orders 
of magnitude greater than that of liquid nitrogen.'^^ Sec- 
ondly, the distortion introduced by the film being em- 
bedded in the substrate, rather than lying upon it, must 
be negligible, since the thickness of the film is far smaller 
than its width or length. Thirdly, the error introduced 
by treating the substrate as semi-infinite must be negligi- 
ble, because substrate and film dimensions typically differ 
by about three orders of magnitude. Fourthly, treating 
substrate and film as being different regions of the same 
infinite piece of homogeneous material means that the 
thermal impedance between the two is zero, whereas the 



accepted value of the actual film-substrate exchange coef- 
ficient is h = 10'^ W/K cmfl^^i^that this approximation 
is acceptable is shown by previous work in which it made 
little difference to the results of numerical calculations.'^ 
Finally, treating substrate and film as being made of the 
same material also means that they have the same dif- 
fusivity; this will be handled by choosing a diffusivity 
coefficient in between that of the real film and the real 
substrate. 



finite 
substrate 



film 

experimental 
setup 




Figure 1. Comparison of the model geometry with the 
schematic experimental setup. In the model, film and sub- 
strate are of the same material so as to be able to use Eq. 
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Splitting the infinite medium 




Figure 2. The model and its mirror image. Heat flow (arrows) 
does not cross the interface between model and image. 

Our next simplification is to assume that the heat- 
ing rate in the film does not depend on position, but 
only on time: Q = Q (t). Since the heating of the films 
we are considering will be caused by electrical current, 
this assumption is equivalent to assuming that the cur- 
rent density is uniform throughout the film. Explicit 
support for this comes from the work of Herrmann et 



al. p'i' l is i who found that in high-Tc tapes current in ex- 
cess of a certain characteristic cutoff (the value of which 
was slightly below the critical current) was indeed dis- 
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tributed homogeneously across the entire superconductor 
cross section. Moreover, direct measurements on a slab of 
Bi-based crystal, using a microarray of Hall probes, show 
that although current is restricted to the lateral regions 
of the slab at relatively low temperature (T < 50 K), 
nonuniformity becomes negligible above about 80 K, a 
temperature relatively near T^^r- Again, a recent stud}*^ — ' 
of YBa2Cu307_5 strips using magneto-optical imaging 
found that even at so low a temperature as 20 K the 
nonuniformity of J is only of the order of 20% when the 
transport current is 90% of its critical value. Finally, a 
uniform distribution of current throughout the cross sec- 
tion is the simplest explanation of the observation that 
critical current is independent of bridge width.'^IHlll 

Given uniform heating and the homogeneity of film and 
substrate (and since no measurements of local tempera- 
tures have yet been made on standard high-Tc bridges), 
we are interested only in the volume-averaged tempera- 
ture of the film. From Eq. ^ , 

T{t) = n + i l^l^l^ ^Ta{r,t\fo,to) dt d% d'r 



( lim M{t — to) = 0) - as in fact was only to be ex- 

fo— ^ — oo 

pected. We exploit this behaviour by making a further 
approximation consisting in the replacement of M as de- 
fined in Eq. ([t]) by Af^., a rectangular function of unit 
height that is non-zero on the interval [0,t], where r is 
the area under M (Fig.|4]): 

m = n+ f ^ M.(i-io) d<o = n+ f ^ dt, 

Jo c Jt_^ c 
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Figure 3. Schematic diagrams of the heating rate Q{to) and 
the transfer function M{t — to) at two times, t and t' > t. 
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where V is the volume of V, A = 2^^ D{t — to) is the 
diffusion length, erf is the error function, and 
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Note that although the precise shape of the transfer 
functiorpSEfil ]\/[ depends on the dimensions and diffusiv- 
ity of the film, its general shape is as shown (reversed) in 
Fig. [3] This means that in the final expression in Eq. ([6]), 
the convolutioiPl of Q{to) with M{t-to) gives consider- 
ably more weight to the immediate past and the present 
moment {\\mM{t — to) — 1) than to the remote past 
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Figure 4. Comparison of the transfer function M{t — to) with 
its approximation, Mr(t — to)- 

We call T the characteristic time of the film. Further- 
more, since Mr clearly overweights the proximal past, 
we slightly correct the integral by taking (5(to) to have 
the value Q{t — t) for all to in [t — r, t], i.e. we replace 
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Mr{t - to) with Ms = Td{t -t -to). Thus, finally (in- 
sofar as the dependence of temperature on time via the 
heating rate), T{t) = TQ{t — t)/c + Th, or 



T{t- 



Q[J,E{J,Tm , ^ 

T \- lb 



(9) 



where the origin of the time-dependence of Q is explicitly 
displayed as the temperature- (and hence time-) depen- 
dence of E{J), the intrinsic current-voltage characteristic 
(CVC) of the film. 

Because of the structure of M as a product of factors 
that each depend on only one spatial dimension (as well 
as on time, through A), r tends to zero whenever w, I 
or d do and the others are held fixed, and tends to the 
product of the other two factors whenever w, I or d tend 
to infinity; the larger w, I and d, the larger is r, i.e. the 
greater the thermal inertia of the film, though t always 
lies in [0,1]. That for a fixed thickness and width/length 
ratio the thermal behaviour of the film does depend on 
width, but increasingly less as width increases, has re- 
cently been confirmed experimentally.^ The dependence 
of T = M {t') dt' on film width and length for a typ- 
ical thickness (0.15 /im) is shown as a contour map in 
Fig.H 



100 



3 




20 



40 60 
w(/um) 



80 



100 



Figure 5. Isolines of the characteristic time r of a typical film 
0.15 fim thick as a function of film width and length. 

Note that if with a fixed current density J, T tends to a 
stable value Too , then Eq. ([9| is asymptotically equivalent 
to Eq. because Q must then also be stable and can 
therefore be taken out from under the integral in Eq. ^ . 

Since our heat source is in this study the Joule effect 
of a constant current, the form of Q[J, E{J, T)] in Eq. (19 



For E{J, r), which is not known (the experimental CVCs 
at constant bath temperature showing the very distortion 
that is investigated in this study), we shall use two em- 
pirical functions with different forms in order to show 
that our results do not depend critically on this aspect 
of the model. The first, varia nts o f which have been used 
in previous work by ourselve^^El] and others ,^221123 is 



EM,T) = Eo{T) [ 



1]' 



(11) 



where Eo{T) = Soi(l - T/T,)™ and Jo(T) = Joi(l 
T/Tc)"\ The second is 



(12) 



where Jo{T) is as in Eq. ([TT]), s = so + si{l-T/T^), and 
Pn (T) is obtained by extrapolation from data for the tem- 
perature dependence of the resistivity of the normal (non- 
superconducting) film. Both have four adjustable param- 
eters _(_Eoi, Joi, mandn in Eq. (11), Jqi, m, Sgandsi in 
Eq. ~ 



Q[J,EiJ,T)]^JEiJ,T) 



(10) 



( |12[ ) and, as in previous work, both are to be fit- 
ted to experimental (J, E) data for the region of small 
E, where little heat is generated and the experimental 
CVC can accordingly be expected to lie close to the in- 
trinsic CVC. We shall call the isotherms corresponding 
to En and Eg, n-isotherms and s-isotherms, respectively. 
It may be noted that whereas Eg tends to the natural 
limit p„ J when J ^ Jq or T — > Tc, this is not so for 

En- 



B. Parameterizing and testing the model 

The films whose thermal behaviour we used to test 
the model and its predictions were two YBa2Cu307_5 
bridges on SrTiOa substrates. One (sample mA) was 
a 50 X 10 X 0.12 /im film with rc=89.8 K, p„(100 
K) = 117^f7cm and p„(300 K)=374uncm; other fea- 
tures have been published elsewhere.'^S The other (sam- 
ple m50a) was longer and wider (500 x 50 x 0.12 ^m), 
with rc=87.1 K p„(100 K) = 190 /Ltllcm, and p„(300 
K) = 490 p,ncmW By way of illustration. Fig. k3^ shows 
the s-isotherms fitted by least squares to experimental 
{J,E) data for sample mSOa. 

The diffusivity D of the superconducting bridge ma- 
terial is 0.05cm^/s, and that of the SrTiOa substrate 
0.18 cm^/s:-^- The intermediate value to be used in the 
model (see the second paragraph of the previous section) 
was informally optimized to afford adequate fit between 
the experimental values of J*(Tft) for sample m50A and 
predicted values that were obtained as follows. 

Although J*{Tb) is defined as the current density at 
which a voltage jump occurs, Eqs. ([9l)-([l2| show that 
this voltage jump will be accompanied by a tempera- 
ture jump. Thus J* (Tf,) may be predicted by using these 
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equations to simulate an experimental plot of tempera- 
ture against current until a discontinuity occurs. Start- 
ing at a bath temperature Tj, at time 0, the temperature 
T(t) attained after applying a current density Jc + SJ 
for time r <C (5i is given by Eqs. ([9|, ([To| and either 



(111 or (12) (where SJ is the current density step used 
in the experiments, typically 0.05 MA/cm^, and St the 
step length, typically 1ms); T(2r) is similarly obtained 
using T(r) as starting temperature; and so on until the 
total time elapsed is St, whereupon the current density 
is increased by SJ, etc. The current density at which 
a sudden temperature jump is observed is identified as 
J*{Tb). 
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Figure 6. Large dots show experimental current-voltage char- 
acteristics (CVCs) of a 50/im-wide YBa2Cu307_i film at bath 
temperatures Tb of (from right to left) 72.1, 73.5, 76.5, 78.0, 
81.2, 82.3 and 84.2 K. Vertical dotted lines indicate volt- 
age jumps at J*{Tb). Continuous curves show s-isotherms 
[Eq. ( 12 1] fitted to the data with E values in the shaded re- 
gion. 

Fig.[7] shows plots of (T — Th) against J obtained in 
this way using s-isotherms, and Tablejl] confirms good fit 
between the predicted values of J* (and those predicted 
similarly using n-isotherms) and the experimental values. 
The achievement of such good fit with such a crude model 
is possibly attributable partly to the fact that for each Tf, 
the value of J*{Tb) lies only about 20% above the largest 
J in the set of {J,E) data to which Eqs. (fTTl) and (12) 



were fitted, i.e. extrapolation was quite limited; partly to 
the temperature and voltage jump occurring only 1-3 K 
above Tb (see Fig.]?]), i.e. far from Tc; and, given these 
circumstances, to the above-noted equivalence of Eqs. ^ 
and ([8]) for r < St and Jc < J < J* ■ 

The value of D affording the above results was 
0.12cm^/s; that this is closer to the difFusivity of the 
subtrate than to that of the superconductor seems rea- 
sonable, since the thermal diffusion length of the sub- 
strate for a time of 1 ms is rather more than 250 /im, 
which is much larger than the film. The same value of 
D also performed well for sample mA. The correspond- 
ing values of r are 0.15 /iS for film mA and 0.55 fis for 



mSOA; as required by the above algorithm, these values 
are both much shorter than the experimental step length, 
1 ms. Note that according to Eq. ([9| the characteristic 
time T is the time taken by the film to return to the 
bath temperature Tf, when heating is stopped. Accord- 
ingly, the use of current pulses lasting just a few tenths 
of a microsecond'^^ ^ and separated by intervals of similar 
length should allow the measurement of current-voltage 
curves that are nearly free from artifactual thermal ef- 
fects. 




J(MA/cm^ 



Figure 7. Excess temperature of sample m50A, relative to 
the bath temperature Tb, as calculated using Eq. ([9| and s- 
isotherms [Eq. ( 12 1] in simulations in which the current den- 
sity J was increased by 0.05MA/cm^ at 1 ms intervals. From 
right to left, Tb = 72.1, 73.5, 76.5, 78.0, 81.2, 82.3 and 84.2 
K. A temperature jump occurs 1-3 K above the bath tem- 
perature, the excess temperature decreasing with increasing 
Tb. 



Table I. Experimentally observed supercritical currents of 
sample m50a at seven bath temperatures, together with the 
values obtained in simulations of experimental runs using 
Eq. (|9| and either n-isotherms [J*; Eq. (Ill or s-isotherms 
j;; Eq. ([T2|]. 
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72.1 


6.27 


6.50 


6.50 


73.5 


5.86 


5.95 


5.95 


76.5 


4.84 


4.80 


4.80 


78.0 


4.04 


4.20 


4.25 


81.2 


2.82 


2.90 


2.95 


82.3 


2.49 


2.40 


2.50 


84.2 


1.61 


1.55 


1.65 
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IV. STABILITY 



Eq. (|9| is a nonlinear autonomous difference equation 
of first orderpS To examine tlie stability of T at constant 
Tb and under a fixed current density J we rewrite this 
equation in the form 



T{t + T)^TQ[T{t)-J] 



(13) 



where 7Q[T(t); J] = r ^S^^^+Tf,, and wc note that un 
der the given conditions Tg (T; J) is convex and monoton- 
ically increasing, and that 7q(T{,; J) > Tf, (since the cur- 
rent must heat the film). The condition for attainment 
of a stable temperature Too is that the graph of Tq{T; J) 
intersect the line T{t + t) = T{t), as may be seen by ex- 
amining the staircase diagrarcPH shown in Fig.js) in which 
vertical arrowed lines indicate real changes in tempera- 
ture during a time increment r, while horizontal lines 
translate the final temperature of one r-interval into the 
starting temperature of the next. Not only does the path 
0^1^-2^3^-- - lead from its starting point (Tf,) to 
the limiting temperature Too (a fixed point oiTqiT; J)), 
but so does the path 1* — 2* — >• 3* — • • • , i.e. any 
fluctuation to a temperature higher than T^o will be re- 
covered from. By contrast, if Tq{T;J) meets the line 
T{t + T) — T{t) tangentially the limiting temperature T* 
is not stable (Fig.joJ) , and if Tq (T; J) lies wholly above 
r(t+r) = T{t) there is no limiting temperature (Fig. 10). 



+ 




Figure 8. Staircase diagram of thermal dynamics under an 
electrical current density J < J* . The inset shows, as a func- 
tion of time, the approach to the stable limiting temperature 
Too from below. 

It may be enlightening to compare the above situations 
with that of a normal conductor, for which E{T, J) = 
p{T)J, where the electrical resistivity p(T) — po + piT. 
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Figure 9. Staircase diagram of thermal dynamics under an 
electrical current density J = J* . Starting at low values, tem- 
perature increases towards the asymptotically limiting value 
T* , but this limit is unstable: starting from higher values, 
temperature increases indefinitely. 
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Figure 10. Staircase diagram of thermal dynamics under an 
electrical current density J > J*. The rate of increase of 
temperature slows down, but finally increases without limit. 



7q(T; J) is in this case a linear function of temperature, 
and whether it intersects the fine T{t + t) — T{t) (i.e. 
whether a stable limiting temperature is attained) there- 
fore depends on the slope of this linear function. For 
copper, for exarrmle, dTg/dT < 1Q~^ for J = 1 MA/cm^ 
and T = 0.14 ps^- so the absence of thermal runaway is 
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ensured. 

Since Tq{T, J) increases with J [by Eq. (10)], increas- 
ing the current density with a given bath temperature 
results in the location of the curve 7q (T; J) progressing 
from that shown in Fig. [8] to that shown in Fig.[lOj and 
examination of the definition of 7q (T; J) shows that the 
greater the thermal inertia of the superconducting film 
(i.e. the greater r), the smaller the current at which in- 
stability sets in. The supercritical current density J* is 
the current density corresponding to Fig.|9j a situation 
that can be characterized by the conditions 



T* = TQiT*;J* 



dT 



iT*;J) = 1 



(14) 



Alternatively, J* can be characterized as the smallest 
current density such that: 



Tq{t-j)>t, vr>Tfc 



(15) 



(and the supercritical temperature T* as the temperature 
at which equality holds) , or as the largest current density 
for which the equation 



rQ{T;j) - T 



(16) 



has a solution. 

To ensure the identification of J* for a given T^, we 
first find the solution Too of Eq. ( 16 1 for a value of J 



just slightly greater than Jc, then for a slightly larger J, 
and so on, until the largest J for which there is a so- 
lution is identified. This procedure is not a simulation 
analogous to those of Section 3, for example, because the 
criterion for increasing J is not the time elapsed but the 
satisfaction of a criterion of convergence that cannot be 
verified experimentally (at least at present). In experi- 
mental practice and simulations, J is not kept fixed in- 
definitely; if J is increased too fast, J* will be missed, 
though the thermal runaway will of course occur (alter- 
natively, if Tq(T;J) approaches T{t + t) — T{t) quite 
closely in the situation of Fig.jlOj the slowing down of 
the increase in temperature may be erroneously taken to 
indicate approach to a non-existent T^). 



For the films considered in Section 3, Fig. 11 shows the 
good fit between the J* {Tb) functions obtained as above 
and the experimental data. Experimental error may rea- 
sonably be regarded as negligible (at least at this repre- 
sentation scale), because the current at which the voltage 
jump takes place is quite well defined, the current-voltage 
curve being locally nearly vertical. That J*(Th) lies at 
higher values for sample mA than for mSOA is expected 
because of its smaller characteristic time and correspond- 
ingly lower 7q(T; J). 



Intuitively, the plots for mA and mSOA in Fig. 11 differ 
because, the narrower the strip of film, the greater the 
proportion of it that is effectively cooled via its edges as 
well as via its bottom surface. It is therefore also to be 
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Figure 11. Predictions of J* for samples mA and m50A ob- 
tained using Eq. Q, s-isotherms [Eq. ( |12[ )] and the conver- 
gence illustrated in Fig. [8] for J values that were successively 
increased until convergence failed. The experimental points 
correspond to the abrupt voltage jumps observed in experi- 
mental CVCs (for sample mSOA, Fig. p|. 



expected that this edge effect will only be significant for 
film strips narrower than a few diffusion lengths. That 
this is so is indeed suggested by Fig. [12] which shows the 
surface J*(T;,,t) calculated using n- isotherms for sam- 
ple mA. For this superconductor and film thickness, only 
films with characteristic times below about 0.05 /iS seem 
likely to be almost free from thermal instability. More 
generally, it appears that other mechanisms that limit 
superconductivity, such as the Larkin-Ovchinnikov elec- 
trodynamic instability, should be studied using films with 
very low r (and/or very short intermittent current pulses, 
as mentioned above) if the effect being studied is not to 
be overwhelmed by the effect of thermal instability. 



< 




Figure 12. J* as a function of bath temperature Tt and the 
characteristic time of the film, r. 
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V. CONCLUDING REMARKS 

This paper has developed a very simple analytical 
model in which strips of superconducting film supported 
by a substrate in a bath thermostatted at temperature 
Tb are represented as rectangular blocks embedded in the 
surface of a semi-infinite medium of the same homoge- 
neous material as the block. Assuming a uniform cur- 
rent density in the block, and given experimental current- 
voltage characteristics that allow parameterization of the 
model, the space-averaged thermal dynamics of the film 
prove to depend on its geometry through a single param- 
eter, its characteristic time r. The quantitative predic- 
tions of the theory as regards the supercritical current 
density J* (Ti,) agree quite satisfactorily with experimen- 
tal observations (Fig. 11). 



Though in this paper the model has been parameter- 
ized for and tested on YBCO films 0.12 /zm thick, its 
derivation involves nothing that prevents its application 
to other high-Tc or low-Tc films, so long as they do 
not require refrigeration with superfluid helium, the ef- 
fective thermal conductivity of which is much greater 
than that of any ceramic substrate-2SJ even for boiling 
films."^ Whether or not the good performance shown 
in Fig.[Tl] also generalizes to other films naturally re- 
mains to be seen. However, the mutual similarity of the 
current-voltage characteristics of all high-Tc supercon- 
ductors suggests that they all become thermally unsta- 
ble under high enough current densities, and that, with 
appropriate parameterization, our model should be able 
to predict this behaviour, thus explaining it in terms of 
thermal instability. 

One of the key assumptions of the model is the homo- 
geneity of the superconducting film, which together with 
the homogeneity of current density guarantees that the 
film and its behaviour are completely described by its ge- 
ometry and current-voltage characteristics. Since YBCO 
films can have various degrees of inhomogeneity (mostly 
in relation to their oxygen content), the successful ap- 
plication of the model to niA and mSOA suggests that the 
inhomogeneity of these samples, if any, was sufficiently 
finely grained as to be negligible. Further investigation 
of this issue probably requires examination of numerous 
individual cases; certainly, it seems safe to suppose that 
inhomogeneities can only accelerate thermal runaway. 

In this paper the intrinsic current-voltage character- 
istics of the film, which are required for calculation of 
its heating rate [Eq. (|lO|)], have been approximated by 
extrapolation from the low-energy region of the experi- 
mental CVCs. Although the experimental CVCs suffer 
from the very inaccuracies, the possible thermal origin 



of which is being investigated, it is assumed that their 
low-energy regions coincide sufficiently closely with the 
required intrinsic CVCs. That the functional form used 
for extrapolation is not excessively critical is supported 
in Table[l]by the agreement between J* and J* and the 
agreement of both with J*^p. Extrapolation will not 
be necessary if ultra-fast current-voltage measurements 
with nanosecond-scale measuring times become available, 
since such measurements may be expected, for the rea- 
sons explained above, to be devoid of thermal distortion. 

The most disconcerting feature of our model is no 
doubt its treating the superconducting film and the sub- 
strate as a single continuum, with a single diffusivity co- 
efficient (the sole free parameter of the model) and no 
acoustic mismatch between film and substrate. How- 
ever, the value of the diffusivity coefficient seems not 
to be excessively critical (a single value worked well at 
all bath temperatures for both the samples considered 
here); while the lack of acoustic mismatch means that 
our results support the significance of thermal instabil- 
ity effects even under the conditions that are least con- 
ducive to such effects, i.e. with optimal thermal coupling 
between film and substrate. Numerous authors have ac- 
knowledged the determinant role of thermal effects when 
thermal impedance is high, but have implicitly or explic- 
itly denied that they are significant when thermal cou- 
pling is good. 

To sum up, in previous studies, finite element 
calculations^ have predicted the occurrence of a ther- 
mally driven transition to the normal conductance state 
at zero applied magnetic field when a homogeneous su- 
perconducting film is subjected to a controlled electrical 
current exceeding a certain "supercritical" value that co- 
incides with experimental observations."^- 

The present work supports those findings analytically. 
Experimental research on the breakdown of superconduc- 
tivity due to other mechanisms (Larkin-Ovchinnikov vor- 
tex instabilities, hot spots, phase-slip centres) must ac- 
cordingly be carried out under conditions that exclude 
the possibility of the heat-driven transition, and vice 
versa. 
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